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Abstract. The bright BL Lac object PKS 2005-489 was 
observed by BeppoSAX on November 1-2, 1998, following 
an active X-ray state detected by RossiXTE. The source, 
detected between 0.1 and 200 keV, was in a very high state 
with a continuum well fitted by a steepening spectrum due 
to synchrotron emission only. Our X-ray spectrum is the 
flattest ever observed for this source. The different X-ray 
spectral slopes and fluxes, as measured by various satel- 
lites, are consistent with relatively little changes of the 
peak frequency of the synchrotron emission, always lo- 
cated below 10 17 Hz. We discuss these results in the frame- 
work of synchrotron self-Compton models. We found that 
for the BeppoSAX observation, the synchrotron peak fre- 
quency is between 10 15 and 2.5 x 10 16 Hz, depending on 
the model assumptions. 

Key words: BL Lacertae objects: general - X-rays: 
galaxies - BL Lacertae objects: individual: PKS 2005-489 



1. Introduction 

The overall spectral energy distribution (SED) of blazars 
(BL Lac objects and violently variable quasars) shows (in 
a v vs vF v representation) two broad emission peaks: the 
lower frequency peak is believed to be produced by syn- 
chrotron emission, while the higher frequency one is prob- 
ably due to the inverse Compton process. The location 
of the synchrotron peak is used to define different classes 
of blazars: HBL (High frequency peaked blazar, with the 
maximum power output in the UV or X-ray frequencies) 
and LBL (Low frequency peaked blazar, peaking in the 
IR or optical bands) (Giommi & Padovani 1994). Thus, 
for LBL sources the X-ray emission is dominated by the 
inverse Compton component. 

Send offprint requests to: G. Tagliaferri 



The continuum emission of blazars is both highly lu- 
minous and rapidly variable and it is probably dominated 
by emission from a jet moving relativistically at small an- 
gles to the line of sight (Blandford & Rees 1978). Strong 
evidence for bulk relativistic motion was first provided 
by multiwavelength studies, later confirmed directly with 
VLBI observations (Vermeulen & Cohen 1994). However, 
single epoch spectra cannot constrain the models of vari- 
ability in relativistic jets (e.g. Kbnigl 1989; Ulrich et al. 
1997). Since blazars emit over the entire electromagnetic 
spectrum, a key for understanding blazar variability is the 
acquisition of several wide band spectra in different lumi- 
nosity states during major flaring episodes. Coupling spec- 
tral and temporal information greatly constrains the jet 
physics, since different models predict different variability 
as a function of wavelength. Important progress in this re- 
spect has been achieved recently for some of the brightest 
and best studied blazars, as PKS 2155-304 (Chiappetti et 
al. 1999; Urry et al. 1997), BL Lac (Bloom et al. 1997), 
3C 279 (Wehrle et al. 1998), Mkn 501 (Pian et al. 1998), 
Mkn 421 (Maraschi et al. 1999; Fossati et al. 2000a,b), 
ON 231 (Tagliaferri et al. 2000). We successfully used the 
BeppoSAX satellite (Boella et al. 1997a) to perform ob- 
servations of blazars that were known to be in a high state 
from observations carried out both in other bands (mainly 
optical and TeV) and in the X-ray band itself. The good 
BeppoSAX sensitivity and spectral resolution over a very 
wide X-ray energy range are ideal to constrain the existing 
models for the X-ray emission. 

PKS 2005-489 (z = 0.071) is a bright BL Lac ob- 
ject that belongs to the 1 Jy radio catalog (Stickel et al. 
1991). Its broad band spectrum peaks in the UV- soft X- 
ray bands making it an HBL type of source. It is one of 
the few extragalactic objects seen in the EUV band (Mar- 
shall et al. 1995) and it has been detected in the 7-ray 
band by EGRET (von Montigny et al. 1995). In the X-ray 
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Fig. 1. The 2-10 keV light curve of PKS 2005-489 during 
the BeppoSAX observation of 1-2 November, 1998. The 
bin time is of 5400 s. 

band it has been observed 5 times with EXOSAT, showing 
large flux variations correlated with changes in the spec- 
trum (harder spectrum with increasing flux) (Giommi et 
al. 1990, Sambruna et al. 1994a, b). These results are also 
confirmed by two ROSAT observations, which also show 
a steep 0.2-2.0 keV spectrum, with photon index of about 
3 (Sambruna et al. 1995). In September 1996, PKS 2005- 
489 has been observed by BeppoSAX, that detected the 
source in a rather bright state, with a flatter 2-10 keV 
photon index of 2.3 (Padovani et al. 1988, 2000). Finally, 
a monitoring campaign has been performed in the period 
October-December 1998 with RossiXTE and the source 
has been detected up to 40 keV (Perlman et al. 1999). 
During this campaign the source underwent a strong flare, 
with a peak to quiescent 2-10 keV flux variation of a fac- 
tor of 8 (30 with respect to the ROSAT fluxes). This event 
essentially lasted for all the RossiXTE campaign, peaking 
on 10 November, thus the source was in a high state for at 
least three months (Perlman et al. 1999). During the week 
of 18-25 October, a flare alert was issued by the all-sky 
monitor (ASM) team (Rcmillard 1998). This triggered our 
ToO observation and PKS 2005-489 was observed with 
BeppoSAX in November 1-2, 1998. 

In this paper we present and discuss the results of this 
BeppoSAX observation together with some simultaneous 
infrared data. 

2. X— ray Observations 

2.1. Observations and Data Reduction 

The BeppoSAX satellite is an Italian mission with the par- 
ticipation of the Netherlands Agency for Aerospace Pro- 
grams (NIVR) and the Space Science Department of the 
European Space Agency (SSD-ESA). It carries on board 



four Narrow Field Instruments (NFI) pointing in the same 
direction and covering a very large energy range from 0.1 
to 300 keV (Boella et al. 1997a). Two of the four instru- 
ments have imaging capability, the Low Energy Concen- 
trator Spectrometer (LECS), sensitive in the range 0.1- 
10 keV (Parmar et al. 1997), and the three Medium En- 
ergy Concentrator Spectrometers (MECS) sensitive in the 
range 1.3-10 keV (Boella et al. 1997b). The LECS and 
three MECS detectors are all Gas Scintillation Propor- 
tional Counters and are at the focus of four identical graz- 
ing incidence X-ray telescopes. The other two are pas- 
sively collimated instruments: the High Pressure Propor- 
tional Counter (HPGSPC), sensitive in the range 4-120 
keV (Manzo et al. 1997) and the Phoswich Detector Sys- 
tem (PDS), sensitive in the range 13-300 keV (Frontera et 
al. 1997). For a full description of the BeppoSAX mission 
see Boella et al. (1997a). 

The data analysis for the LECS and MECS instru- 
ments was based on the linearized, cleaned event files ob- 
tained from the online archive (Giommi & Fiore 1998). 
Light curves and spectra were accumulated with the 
FTOOLS package (v. 4.0), using an extraction region of 
8.5 and 4arcmin radius for the LECS and MECS, re- 
spectively. At low energies the LECS has a broader Point 
Spread Function (PSF) than the MECS, while above 2 keV 
the PSFs are similar. The adopted regions provide more 
than 90% of the source counts at all energies both for the 
LECS and MECS. 

As expected, PKS 2005-489 was detected in a high 
state in the X-ray band. However, our observation oc- 
curred about one week before the source reached its max- 
imum as monitored by RossiXTF* (Perlman et al. 1999). 
Some flux variability of the order of 10% on a timescale 
of about 10 hours was detected during our observations in 
the 2-10 kev energy band (see Fig. 1). No similar flux vari- 
ability was detected below 1 keV. Given that the amount 
of variability is quite small, we extracted only one spec- 
trum for the full observation. This corresponds to a total 
LECS exposure time of 20395 s and MECS exposure time 
of 52478 s. The background files for the spectral analy- 
sis were accumulated from the long blank field exposures 
available from the SDC public ftp site (see Fiore et al. 
1999, Parmar et al. 1999). 

The HPGSPC and PDS were operated in the custom- 
ary collimator rocking mode, with each collimator point- 
ing alternatively at the source and at the background for 
96 s. The HPSGPC has a single collimator, thus the whole 
instrument is looking at the target for half of the time 
and to an offset direction for the other half. The total ex- 
posure time in our case was 23258 s. The PDS has two 
collimators, thus one is always looking at the source and 
the other one at the background. The total PDS exposure 
time, for the two collimators, was of about 46.8 ks. The 
background subtracted HPGSPC and PDS spectra were 
obtained from the standard pipeline analysis carried out 
at the BeppoSAX Science Data Center. 
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2.2. Spectral Analysis 

For the spectral analysis, the LECS data have been consid- 
ered only in the range 0.1-4 keV, due to still some unsolved 
calibration problems at higher energies. To fit the LECS, 
MECS, HPGSPC and PDS spectra together, one has to 
introduce constant rescaling factors to account for uncer- 
tainties in the inter-calibration of the instruments. The 
acceptable values for these constants are in the range 0.7- 
1.0 for the LECS and in the range 0.77-0.93 for the PDS, 
with respect to the MECS. The HPGSPC constant is nor- 
mally in the range 0.9-1 (Fiore et al. 1999). The spectral 
analysis was performed with the XSPEC 10.0 package. 

In the fitting procedure we first considered a single 
power law model plus absorption with the column den- 
sity fixed at the Galactic value Ah = 4.2 x 10 20 cm~ 2 . 
This model does not give a good fit to the data yielding 
a reduced \r — 2.0 (214 degree of freedom). The fit im- 
proved by allowing A^h to vary (see Tab. 1). The power 
law fit gives a good representation of the 0.1-200 keV X- 
ray spectrum, with a relative steep photon spectral in- 
dex a = 2.18. This index is slightly flatter than the one 
measured by RossiXTE (a = 2.37 ± 0.01 in the 2-10 keV 
energy band) in the same days (1-2 November, Perlman 
et al 1999). Actually, the spectral index determined with 
BeppoSAX is flatter than all spectral indices obtained with 
RossiXTE in the period October December, 1998. The 
flattest RossiXTE spectral index is a = 2.32 ±0.015) mea- 
sured on November 6, 1998, a few days before the top of 
the X-ray flare (Perlman et al 1999). Moreover, the very 
high statistics of the LECS and MECS spectra show that 
from a statistical point of view the power law fit is not ac- 
ceptable and that a more complex model is required. This 
is also reflected by the fact that the interstellar absorption, 
Ah, determined in the best fit procedure is higher, at 90% 
confidence, than the Galactic value. This implies either 
the existence of an intrinsic absorption in PKS 2005-489 
which seems unlikely, or that we are using a wrong model. 

We therefore considered a broken power law model 
which gave a much better fit, even if we keep the fixed 
at the Galactic value. The difference between the two spec- 
tral indices, although significant, is small, with the break 
at ~ 2 keV. The spectral index above the break is the same 
as the single power law, while before the break it is some- 
what flatter (see Tab. 1 and Fig. 2). Of course, this break 
at 2 keV could not have been detected by RossiXTE. In 
this fit the LECS, HPGSPC and PDS intercalibration con- 
stant factors have also better values (see notes in Tab. 1). 
Clearly we are detecting a downward curved spectrum, 
probably the declining part of the synchrotron emission. 
Thus, we considered a more physical intrinsically curved 
spectral model developed by Fossati et al (2000b). This 
model is an analytic expression that assumes a continu- 
ously changing slope between two asymptotic limits. It is 
expressed in a form that provides the following parame- 
ters (£1, ai, E 2 , Q.2, Eb, /), where E\ and E2 are the two 
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Fig. 2. LECS+MECS+HPGSPC+PDS spectrum ob- 
served on 1-2 November 1998, the best fit model is a bro- 
ken power law model. The N-r was fixed to the Galactic 
value. 

"pivoting energies", while Eb and / determine the scale 
length of the curvature. In order to fit this model to the 
data we have to fix one parameter for each of the two 
pairs (Ei, a\) and (E 2l a 2 ). If we fix the energy we will 
find the spectral index at that energy, while if we fix the 
spectral index we will find the energies where it will oc- 
cur. The latter possibility is very important, because it 
allows us to find the peak of the synchrotron component, 
if this is inside our energy band, and to estimate its un- 
certainty. This can be obtained by setting one spectral 
index, for instance a\ = 1, so that the best-fit value of E\ 
will give Ep ea k. For more details we refer to the Fossati 
et al. (2000b) paper. In order to determine the values of 
the really interesting parameters and their confidence in- 
tervals, in the best-fit procedure we kept the absorption 
column density fixed to the Galactic value and / fixed 
to the best-fit value / = 4.5. If we fix the two energies 
at Ei = 0.1 and E 2 = 100 keV, we find that the re- 
spective photon spectral indices are ct\ = 1.987_o!o2 and 
a 2 = 2.23 ± 0.02. We then fixed the photon energy index 
a 1 = 2 in order to find the synchrotron peak and indeed 
we find that E pea k = 1.157_o.'i5 keV. This confirms the 
presence of a curved spectrum with a synchrotron peak 
below <~ 1.5 keV. It is in any case remarkable that this 
curvature from 0.1 up to 100 keV is rather small with re- 
spect to other sources. All quoted errors, for the curved 
model, are calculated for A% 2 = 2.7 that corresponds to 
a 90% confidence range for one parameter of interest. 

3. Infrared Observations 

The JHKL observations were made using the Mklll in- 
frared photometer attached to the SAAO 1.9m telescope 
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Table 1. Fit results for a power law and broken power law model 



model 


I\ b or N H 


r 2 


break 
energy c 


X 2 (d.o.f.) 


r [0.1-2 keV] 

erg cm~ 2 s _1 


-F[2-10 keV] 

erg cm -2 s _1 


-^[10-100 keV] 

erg cm~ 2 s _1 


power law 
broken p. 1. 


(4.9;S:iS) * io 20 

2 02" 04 

z.o^ +0 03 


9 i O-0.02 
z - lo +0.02 
9 91 -0.02 
z - zi +0.02 


1 or-0.35 
1.8t> +0 45 


1.4 (213) 
0.90 (212) 


4.0 x 10~ 10 


1.8 x 10~ 10 


1.7 x 10~ 10 



a we had also a constant to allow for uncertainties in the intercalibration of the instruments (see sec. 2.2); the best fit values 
are 1) PL: LECS constant= 0.75, MECS constant= 1. frozen, HPGSPC constant^ 0.94, PDS constant= 0.72; 2) BPL: LECS 
constant^ 0.78, MECS constant= 1. frozen, HPGSPC constant= 0.97, PDS constant^ 0.77; b photon spectral index; c energy 
values in keV; d corrected for the absorption. 



Table 2. Infrared Observations 



date 


Jul. Date 


J a 


H 


K 


L 


27-07-97 


2450657 


10.97 


10.23 


9.55 


8.52±0.03 




2450721 


11.17 


10.41 


9.75 


8.80±0.08 




2450738 


11.23 


10.51 


9.87 


8.91±0.05 




2450742 


11.26 


10.49 


9.84 


8.70±0.07 


30-10-98 


2451117.3 


11.16 


10.46 


9.80 


8.69±0.04 



a the errors for the J, H, K magnitudes are < 0.03. 



at Sutherland, in July 1997 and in October 1998, the latter 
was one day before the start of the BeppoSAX observation. 
The chopper throw was 30 arcsec in a N-S direction and 
the aperture was 12 arcsec diameter. This aperture contain 
a large fraction of the galaxy whose H-flux is comparable 
(<~ 20% lower) to the H-flux of the nucleus (see Kotilainen 
et al. 1998). The values reported in Tab. 2 and Fig. 4 are 
not corrected for the galaxy contribution. However, on the 
scale of Fig. 4 this correction is negligible and smaller than 
the size of the infrared data symbols. The calibration in 
terms of log flux in W m -2 Hz -1 for a zeroth mag star on 
this system at JHKL arc -22.82, -23.01, -23.21 and -23.55 
respectively. 

As it is usually found for HBL, there is very little 
change in the slope of the infrared spectrum as the lu- 
minosity changes (see Tab. 2), i.e. the infrared colour of 
the blazar changes little or not at all. This is due to the 
fact that the synchrotron peak is at much higher ener- 
gies. Changes in K of up to four mags have been seen in 
some blazars with little change in their infrared colours 
(e.g. Fan 1999; Hagen-Thorn et al. 1994). In the partic- 
ular case of PKS 2005-489, the K mag changed only by 
0.32 mag during the observations, which is relatively little 
for this class of objects. 

4. Discussion 

4.1. SED 

We observed PKS 2005-489 while it was in a very high 
state, aiming at detecting a change in its spectral energy 
distribution (SED) with respect to previous observations. 
Very significant changes of the overall SED were in fact 



detected in other HBL objects, such as MKN 501 (Pian et 
al. 1998) and 1ES 2344+514 (Giommi et al 2000): in these 
objects, during flares, the synchrotron peak frequency in- 
creased by 2-3 orders of magnitudes with respect to qui- 
escent states. These observations provide new conditions 
that the various models have to reproduce, exploring new 
range of values for the physical parameters and allowing 
us to investigate the accelerating mechanisms inside the 
jets. 

In the case of PKS 2005-489, the X-ray data can be 
fitted equally well by a broken power law and by a continu- 
ously steepening model. Although the spectral parameters 
are not very different, the two models predict a quite dif- 
ferent frequency for the peak of the synchrotron emission, 
as illustrated in the two panels of Fig. 4. In the case of a 
broken power law in fact, it can be possible to interpret 
the entire far-IR to hard X-ray SED as due to a single 
synchrotron component, peaking at <~ 10 15 Hz. If, instead, 
the X-ray data are fitted by the continuously steepening 
model, then the synchrotron peak is at ~ 0.1 keV, and the 
IR flux must be produced by another jet component. This 
uncertainty in locating the synchrotron peak frequency 
originated the two possible sets of parameters of the ho- 
mogeneous synchrotron self Compton model by which we 
can interpret the overall SED, discussed below. 

In any case, our data show that the source was in a very 
high state and that the 0.1-200 keV spectrum is well fitted 
by a spectrum due to synchrotron emission only. In Fig. 3 
we plot an enlargement of the SED in the X-ray spectral 
region. It can be seen that PKS 2005-489 underwent dra- 
matic flux changes, but, at least above ~1 keV, the spec- 
tral shape was always steeper than unity, constraining the 
peak of the synchrotron emission to be at lower energies. 
Our X-ray spectrum is the flattest ever observed, even 
considering the many observations by RossiXTE during 
the same period of our 1998 BeppoSAX observation (Perl- 
man et al. 1999). 

4-2. Limits on magnetic field and particle energies 

The spectral properties of our X-ray data indicate that the 
dominating radiation process is synchrotron. Therefore 
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Fig. 3. The X-ray energy spectra of PKS 2005-489 as 
observed at different epochs by different satellites. Source 
of data: Sambruna et al 1994b (EXOSAT), Comastri et 
al. 1997 (ROSAT), Perlman et al. 1999 (i?oss?XTE). 



the emission is due to the most energetic electrons, suf- 
fering rapid radiation energy losses. During our observa- 
tion the flux variability was modest (10% on a timescale of 
about 10 hours), but EXOSAT, in two occasions, observed 
flux changes of a factor of about 4 in 4-5 hours (Giommi et 
al. 1990). We can then require that the emitting electrons 
cool radiatively (by synchrotron and inverse Compton) on 
these timescales, implying a lower limit on the magnetic 
field. Assuming equal radiation and magnetic energy den- 
sities, using i V ar = 4 hours and v x = 1 x 10 17 Hz, we 
derive B > O.ZS' 1 / 3 Gauss and j x < 3 x lO 5 ^ 1 / 3 . Here 
S = [r — Vr 2 — 1 cos is the Doppler beaming factor, 
where 9 is the viewing angle. 



Given the uncertainty in the location of the syn- 
chrotron peak, we have considered two models, shown in 
the two panels of Fig. 4. The top panel shows the X-ray 
data fitted by a broken power law model, and in this case 
it is possible that the simultaneous IR flux is produced in 
the same region producing the high energy emission. We 
have then considered a source with size R — 1.5 x 10 16 
cm and a beaming factor S = 18, resulting in a minimum 
observed variability timescale of ^6 hours, which is typ- 
ical for this class of sources, and close to the minimum 
variability timescale observed in this source by EXOSAT 
(Giommi et al. 1990). With these values of the size and 
beaming factor, the magnetic field is constrained to be of 
the order of B — 1 Gauss, in order to have nearly equal 
synchrotron and self-Compton luminosities, as observed. 
The intrinsic power (assumed to be continuously injected 
throughout the source in the form of relativistic electrons) 
is required to be Li n j = 5.5 x 10 41 erg s _1 . The relativistic 
electrons are assumed to be injected with a power law en- 
ergy distribution of slope s = 2.2 between jmm = 3 x 10 3 
and j max — 10 6 . The value s — 2.2 is required to match 
the slope of the X-ray spectrum. The resulting SSC spec- 
trum is shown as a solid line in the top panel of Fig. 4. 
It can be seen that this model also agrees with the (ad- 
mittedly uncertain) flux in the EGRET band, and well 
describes the entire SED (except for the mm-radio part, 
which should come from larger jet regions, less affected by 
synchrotron self-absorption at these frequencies). Being 
produced by the same region, we expect some variability 
correlations between the IR and X-ray flux. In particular, 
if the radiative cooling time at both frequencies is shorter 
than the light crossing time, then it is the latter timescale 
that is observed in both energy bands. With the adopted 
parameters, this is indeed the case: at the observed fre- 
quency v = 3 x 10 14 Hz, the intrinsic cooling time is half 
the crossing time. 



4-3. Homogeneous SSC model 

We have fitted the SED of PKS 2005-489 with a pure, 
homogeneous and one zone SSC model, assuming no con- 
tribution from the photons produced externally to the jet 
as seeds for the Compton scattering process. This class of 
SSC models is completely constrained once the peak ener- 
gies and their fluxes are known. Unfortunately, in our case 
only the X-ray and IR data are simultaneous, and the lo- 
cation of both peaks (particularly the Compton one), is 
somewhat uncertain. However, the X-ray data alone con- 
strain the synchrotron peak to be at or below the low 
energy end of the BeppoSAX range (i.e. < 0.1 keV). For 
the high energy part of the SED we have considered the 
previous observation by EGRET, resulting in a marginal 
detection, see Fig. 4, (Lin et al. 1997; von Montigny et 
al. 1995) and the TeV upper limit given by CANGAROO 
(Roberts et al. 1999). These data constrain the Compton 
peak to be close to 10 25 Hz. 



The bottom panel of Fig. 4 shows the X-ray data fitted 
by the continuously steepening model. In this case the IR 
flux cannot be produced by the same region emitting the 
X-rays, resulting in a significantly larger frequency of the 
synchrotron peak. The adopted parameters of the SSC 
model are similar to the previous case (see the caption 
of Fig. 4). The only significant difference is the value of 
Imin = 1-5 x 10 4 , now much larger. 

In this case we predict /expect that the IR flux has a 
different variability behavior than the X-ray flux, since it 
would probably be produced in a larger region of the jet, 
with a smaller value of the magnetic field, making the cool- 
ing time of the IR emitting electrons longer. Furthermore, 
if the variability is caused by a perturbation of a steady jet 
propagating from the inner to the outer regions, we expect 
some delay between the X-ray and the IR flux variations, 
giving information about the relative distance of the two 
regions. 
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Fig. 4. The SED of PKS 2005-489. Note that only the BeppoSAX and the IR data (filled symbols) are simultaneous. 
The other data have been taken from the Nasa Extragalactic Database (NED). In the top panel we used for the X-ray 
spectrum the broken power law best-fit model. In this case the SSC model fits also the IR data. In the bottom panel 
we used for the X-ray spectrum the curved best-fit model (see text). In this case the IR emission is not due to the 
same SSC component responsible for the X-ray emission. For both SSC models we have assumed no external photons 
contributing to the inverse Compton process. The input parameters for the top (bottom) models are: L' in j = 5.5 x 10 41 
(7.8 x 10 41 ) erg s~\ B = 1 (1) Gauss, R = 1.5 x 10 16 (10 16 ) cm, 7mm = 3 x 10 3 (1.5 x 10 4 ), lmax = 10 6 (1.25 x 10 6 ), 
s = 2.2 (2.3), 8 = 18 (15). 



